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Abstract 

The final-state interaction in multichannel decay processes is systematically studied in the 
hadronic picture with application to B decay in mind. Since the final-state interaction is intrinsi- 
cally interwoven with the decay interaction in this case, no simple phase theorem like "Watson's 
theorem" holds for experimentally observed final states. We first solve exactly the two-channel 
problem as a toy model in order to clarify the issues. The constraints of the two-channel approxi- 
mation turns out to be too stringent for most B decay modes, but realistic multichannel problems 
are too complex for useful quantitative analysis at present. To alleviate the stringent constraints 
of the two-body problem and to cope with complexity beyond it, we introduce a method of ap- 
proximation that is applicable to the case where one prominent inelastic channel dominates over 
all others. We illustrate this approximation method with the amplitude of the decay B — > Kn fed 
by the intermediate states of a charmed- meson pair. Even with our approximation we need more 
accurate information of strong interactions than we have now. Nonetheless we are able to obtain 
some insight in the issue and draw useful conclusions on general features on the strong phases. 
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I. INTRODUCTION 



The well-known phase theorem holds for the final-state interaction (FSI) of decay 
processes when the final state consists of a single eigenstate of scattering. While no simple 
nontrivial extension is known in the case of multichannel final states, some calculations were 
made in the past with unjustified extension of the single- channel phase theoremi2]. A two- 
channel problem was studied with a certain class of S-matrix and the correct observation 
was made that inelastic channels are the main source of strong phases in many B decay 
modes[§]. However, it is not easy to obtain quantitatively rehable results from the two- 
channel model. Taking the large limit of open channels, the statistical model 0, was 
proposed as an alternative approach. Quantitatively, however, it is short of predictive power 
since it does not ask for detailed knowledge of strong interaction. When one approaches the 
problem in the quark-gluon picture, one faces inability or large uncertainty in computing 
contributions of the soft collinear constituents numerically. It is fair to say that at present 
we are far from successful computation of a FSI phase in multichannel decay processes. 

In this paper we first study the two-channel problem in detail. The problem is solvable in 
a reasonably compact form without approximation or assumption if relevant information is 
available about strong interaction physics and decay branching fractions. The general solu- 
tion to the two-channel toy-model shows how the elastic scattering phases and the channel 
coupling contribute to the total FSI phase. It concludes in agreement with Donoghue et al^ 
that if a large strong phase emerges in the B decay into two light-mesons, its major source 
is coupling to decay channels that have large branching fractions.^ Although it points to 
the source of problems in strong phases, the simple two-channel model is inapplicable to B 
decay. We proceed to the case of more than two decay channels. Even the three-channel 
problem is mathematically too complicated for solving in a compact form. On the physics 
side our knowledge of strong interaction at total energy 5 GeV (~ m^) is not good enough 
to carry through the analysis with precision. To cope with formidable complexity of the 
problem, we introduce an approximation method that works in the case that, aside from 
the channel of our interest, one inelastic channel dominates over all others. This is different 
from an approximate two-channel problem. It can happen, for instance, to the two-body 
light-hadron decays of B meson when they couple to the charmed-meson pair states. We 
apply our approximation method to the decay B — > Kit and make semiquantitative analysis 
with knowledge of hadron physics currently available to us. 



II. FRAMEWORK AND INPUT 

Two basic ingredients in discussion of FSI are unitarity and time-reversal. In the Standard 
Model the decay interaction is sum of effective local operators Oaia = 1,2,3, ■ ■ ■), each of 
which has a CP-violating (and therefore T-violating) phase factored out as 

H,,,te''- +h.c., {TH,^tT-' =H,^t)- (1) 

The T-violating "weak phase" 6^ arises from the CKM elements. In computing decay 
amplitudes we work separately on the T-invariant part Hint of the decay interaction. In the 



This was suggested to the author by Wolfenstein in many occasions over years. 
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case that T-violation is more general, we can break up the weak interactions into T-even 
and T-odd parts as H■^^ = \{Hu, ± TH^T'^) so that 

TH^T-' = ±H^. (2) 

Then both h\^^ and iH\^^ are T-even since i ^ i* = —i under time reversal.^ We should 
compute the decay amplitudes for H^^^ and iH^^^ separately and take a suitable sum of them 
at the end. Therefore it is sufficient to consider in general only T-even weak interactions. 
No matter what method one may use, computation of FSI phases must always be made 
separately for different decay operators. Because two decay operators generate two different 
FSI phases for the same decay process even if the net quantum numbers of operators are 
identical. We consider in this paper only final-state interaction of strong interaction though 
it is in principle easy to include electromagnetic FSI. We shall refer to FSI phases also as 
strong phases in this paper. 

With T-invariance, the strong ^-matrix operator obeys 

TST-^ = Si (3) 

We can always choose phases of states such that the T-invariant S'-matrix elements [Skj = 
{k\S\j) = are not only unitary but symmetric; 

Sjk = Skj (4) 

since — > and {k°^'^\ \k^^) under time reversal. It is emphasized here that 

the requirement of Eq. (jlj) fixes the phases of states except for the overall sign of ±1.^ 
Specifically for the eigenchannels of the S'-matrix |a™) and it holds by definition that 

{b\S\a) = {b'"''\a''') = 6bae^'^\ (5) 

where Sa is the eigenphase shift. In the case of decay matrix elements, the initial state is a 
one-particle state that is stable with respect to strong interaction. Since the initial decaying 
state is an asymptotic state with respect to strong interaction, there is no distinction between 
"in" and "out" states. For B decay 

\B) ^ {Bl (6) 

where we choose the i? meson at rest. Eq. removes an arbitrary unphysical phase from 
the state \B) too. 

A simple relation results from time reversal of the decay matrix element {a°^^\Hint\B) 
when the final state (a°"*| is an eigenstate of S. T-invariance of Hint leads to 

(a°"*|i/,„i|i?) = {B\H,nt\an. (7) 



^ If we make this breakup for the Standard Model interaction, we would get = cos S^Hint and 

iH^~^ — — sinJw-ffmt for the first term Hy^ — Hinte^^'" of Eq. ([T]). 

If one multiplies the states with some phases as |j) — > e*"|j') and \k) e^^\k), Sjk and Skj would acquire 
phases of opposite signs e'^^'^^~"^ so that the equality Sjk = Skj would break down. 
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Inserting the completeness relation J2b = I next to Hint in the right-hand side, 

we obtain from Eq. ([7j) with Eq. ([5]) 

(a°"*|i/,„t|i?) = ^(5|i/,„i|6°"*)(6°"*|a™) = e''''^{a°^'\H,nt\By. (8) 

b 

With the decay amplitude {a°'^^\Hint\B) into the eigenchannel a denoted by Aa, this relation 
reads 

A, = e^^'^Al, (9) 

namely, Aa = ie**^" I | ■ This is the well-known phase theorem usually referred to as Wat- 
son's theoremil,]. It is a powerful theorem when the final state is an eigenstate of S. It has 
been an important tool of analysis in hyperon decay and K ^ irn decay. 

But this relation is of little use when rescattering has inelasticity. We mean by inelastic- 
ity that observable final states are linear combinations of eigenstates whose weights are not 
simply determined by the Clebsch-Gordan coefficients of isospin symmetry or SU(3) coeffi- 
cients. In such cases the phase of the elastic scattering amplitude has little to nothing to do 
with the FSI phase of the corresponding decay amplitude, as we shall see it in a moment. 

Usefulness of the phase theorem is thus limited to the decay of low-mass particles where 
rescattering is purely elastic up to isospin structure. If the K meson mass were sufficiently 
above 1 GeV, for instance, vrvr of definite isospin would no longer be an eigenstate of S'-matrix 
even approximately. The state pp and uuj would enter an S'-matrix eigenchannel with tttt 
and composition of such an eigenstate depends on low-energy dynamics of the transition 
among vrvr, pp, and uu. In the case of B decay an experimentally observed final state is a 
linear combination of many different S'-matrix eigenstates so that the net FSI phase results 
from the eigenphases weighted with the decay amplitudes of S — >■ eigenstate. Take, for 
instance, the vrvr final state in J = of 5° decay. The state |7r7r)7=o (in s-wave) is far 
from being an S- matrix eigenstate at energy ttib- If we want to use the phase theorem Eq. 
iQ, we must expand Ivrvr) in the strong S'-matrix eigenstates at m^. However, we have 
little knowledge of these eigenstates since their composition depends sensitively on strong 
interaction at long and intermediate distances. Experimentally the two-body channels do 
not account for all final states in B decay. Three and four particle final states of the same 

jPC 

may be significant. Unless there is a good reason to believe that channel coupling is 
negligible among these final states, the strong S'-matrix eigenstates at are made of many 
different particle states (vrvr, pp, KK, ttttKK, • • •). Therefore, if we expand the tttt state at 
total energy m^, it is a linear combination of many different eigenstates of strong S'-matrix. 
We would have to know these expansion coefficients in order to determine the FSI phase of 
TTTT final state with the phase theorem. 

Let us formulate what we have described above. When an observable final state \i°'^'^) is 
not S-matrix eigenstate, we expand it in the eigenstates of S-matrix |a°"*) as^ 

|^-t) = ^0,Ja-*). (10) 

a 

We choose that the S'-matrix is symmetric {Sij = Sji) in the basis of the observable states, 
(cf Eq./ (jl])) Then the expansion coefficients Oia are real, that is, the transformation matrix 



^ We represent the S'-matrix eignestates by \a), \b), - ■ ■ and the observed particle states by |fc), • • •. 
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O in Eq. (fTOi) is an orthogonal matrix. (See Appendix if proof is needed.) With this 
expansion and Eq. Qj, 

A = E 0^aAa = E O.^e'^'" (11) 
a a 

We are able to write the right-hand side of Eq. (ITT!) in terms of observable decay amplitudes 
Aj s as 

= T.S^,^P (12) 

j 

where Sij = OiaC^^^'^Oja has been used. This is the fundamental relation in discussion of FSI. 
The physical picture is simple and clear. (Fig.l) The input is unitarity and time- reversal 
of the S'-matrix aside from choice of unphysical phases of states. Dynamical information of 
strong interaction is fed through the eigenphases 5a and the orthogonal mixing matrix O. 
In addition, we must provide relative magnitude of the amplitudes Ai as independent pieces 
of input from weak interaction. Consequently the phase of Ai depends not only on strong 
interaction but also on weak interaction. It is clear here that the FSI phase of a decay 
amplitude has virtually nothing to do with the phase of the elastic scattering amplitude 
a j{s) of J = at 5 GeV. Although it looks almost futile to go any further, the purpose of 
this paper is to extract something useful for B decay out of Eq. f|T2|) . 




FIG. 1: The final-state interaction relation in diagrams. Eq. (12). 



III. DYNAMICAL INPUT 

We focus on the FSI phases of the two-body decay modes of the B meson. The phases 
of three-body decay amplitudes depend on the sub-energies of three particles. It is only the 
phases integrated over the sub-energies with the total energy fixed to that enter Eq. 
( fT2|) . We need the S'-matrix elements Sij from experiment. To be concrete, let us consider 
the elastic 7c~^tt~ scattering amplitude as an example. The argument below is identical for 
other two-body channels of light mesons. In B decay the relevant partial-wave channel is 
7r~^n~ in s-wave {J^ = O"*") with isospin / = and 1 = 2. 

Although the high-energy vrvr scattering cannot be directly measured in experiment, we 
can make a reasonable estimate about elastic vrvr scattering since the center-of-momentum 
energy 5 GeV is in the high-energy asymptotic region well above the vrvr resonances.^ The 



^ The excited charmonia exist at mass not far below 5 GeV but their coupUng to Ught hadrons is suppressed 
by QCD. CoupUng of tttt to the open charm channels^] is one of interesting subjects of our study later in 
this paper. 
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FIG. 2: The Regge exchange in elastic vrvr scattering. 

two-body light-hadron scattering in the high energy asymptotic region was studied theo- 
retically and experimentally in the 1960 's. The Regge theory describes this physics well. 
The properties of the Regge trajectories can be deduced from meson-nucleon and nucleon- 
nucleon scattering even though we have no meson- meson scattering experiment. We give a 
very brief review [8|] of our knowledge in this area of the past time since it is an important 
input in our study of FSI. 

First of all, the Pomeron exchange dominates in high-energy elastic scattering. (Fig 2) 
The Pomeron may include a cut and can be more a complicated singularity than a simple pole 
in J-plane unlike the non-Pomeron trajectories such as p and f2- At the level of numerical 
accuracy of our discussion, however, we treat the Pomeron as a simple pole at J = a{t) 
with the intercept a(0) = 1 and a vanishingly small slope. This entails factorization of 
the J-plane residue into product of two vertices. Since isospin is zero for the Pomeron, it 
contributes equally to the / = and 1 = 2 states of the crossed channels tttt. With these 
properties of the Pomeron we obtain necessary pieces of information on tttt scattering from 
TTp and pp scattering at high energies. 

The invariant amplitude for the asymptotic elastic scattering is parametrized with the 
Regge parameters of the Pomeron in the form 



sin Trap (t) ^Sq^ 



(13) 



where Pp^'it) = {Pp{t)y and ap(t) ~ 1 + a'p{0)t + ■ ■ ■ are the factorized residue and the 
trajectory of the Pomeron, respectively, in terms of the invariant momentum transfer t.^ 
We can fix /3p'^(0) by the optical theorem oVot = ImA(s, 0)/s and the factorization relation 



'tot 



Wtot) fJl 



pp 

ot 



22mb at s/s 



TTTT 

tot 



2 7rp 
^'^tot 



rriB 1? 

, (2\2 PP[ 

- UJ '-'totx 



^The value 22mh is in line with the empirical 
The width of the forward peak is more or 
We approximate the Pomeron slope 
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quark counting rule, cx- 

less universal to all elastic hadron scatterings [111, 
to zero (ap(0) ^ 0) and fit the forward elastic peak with the standard exponential form 
exp(t/to)- The diffraction peak parameter is fixed by experiment to to — (0.22 ~ 0.29)GeV^ 



by the elasticity (Jeij (Jtot at i/i = mslol, |13|, |l4]. This value of to reproduces the t dependence 
of daei/dt that falls by r oug hly three orders of magnitude from t = to t ~ — 1 GeV^ for 



pp and Tip scattering [ll|, |12|. It is justified by factorization to choose the vrvr forward peak 



parameter equal to that of pp and vrp. The vrvr invariant amplitude at high energies is thus 



^ The possible small log^ s rise of atot{s) does not set in at y/s ~5 GeV. We may safely ignore it in our 
numerical work. We have included in (3p'^{t) the i-dependence of some other kinematical factors such as 
2ap{t) + 1. The value of sq is at our choice, normally chosen to be 1 GeV^. 
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set to 

~ 22 mb X 2 s ei, (14) 

where we choose to = (0.253±0.033)GeV^. The uncertainty in is primarily due to whether 
one estimates it with aei/cTtot of vr^p or pp and pp. 

The partial-wave amplitudes a/(s) can be projected out of A'^'^{s,t). The result for the 
s-wave is: 

aj^(s) = (0.282 ± 0.037)z, (at s = m|), (15) 
which leads to the s-wave S'- matrix with 2\p^^\/^/s ^ 1, 

S'o""(s) = l + 2mr(s). (16) 

Hereafter we shall often parametrize strength of elastic scattering by e as 

So{s) = 1 - e, (17) 

With Eq. ( |T5i) . the Pomeron contribution to the S'-matrix of tttt scattering at is 

S^"" - 1 - (0.56 ± 0.07). (18) 

The partial- wave amplitudes ai{s) extracted from the flat Pomeron amplitude are purely 
imaginary for all /. The amplitude ao(s) approaches asymptotically the imaginary axis 
below the center of the Argand diagram (shown schematically for / = in Fig. 3). If one 
described the high inelasticity of high-energy vrvr scattering by an absorptive black sphere 
potential, one would have Sq'"{s) {i.e., ao(s) 0.5i). In this limit ad |crtot for 
all Vs by shadow scattering effect, which is in disagreement with experiment. Although the 
numerical value in the right-hand side of Eq. (JTSll has been extracted for the vrvr channel, it 
is much the same for other two-meson channels. With a help of the Kp cross section cXtof 
we obtain 

S^^ - 1-0.51, (ttK) 

5"^^- 1-0.45, (KK), (19) 

where uncertainties are comparable to ±0.07 quoted for vrvr in Eq. (ITSl) or a little larger. 

The values in Eqs. (fTSjl and (fT9!) are the Pomeron contribution alone. The nonleading 
Regge exchanges generate a small imaginary part for Sq. The relevant trajectories are the 
p and /2 in the case of tttt. Their contributions can be estimated with a few additional 



theoretical inputs^ from the cross section difference a^-p — cr^+p — 1.6m6 at ^/s = mp|l4 
Within the uncertainty due to t-dependence of the residue Pit), their contributions to ImS'o 
are at the level of 0.05i for vrvr. Some may wonder about validity of extracting the s-wave 
amplitude from the forward peak region alone. The s-wave amplitude has a flat angular 
dependence so that the contribution of ao(s) extends equally to all directions (Po(cos 9) = 1). 
On the other hand experiment shows that the forward peak falls off by more than three 
orders of magnitude and there is no sign of the s-wave contribution at large angles. But 
this is no surprise. The s-wave amplitude at large angles is canceled by the partial-wave 
amplitudes of up to Z = 0{^/s) which are rapidly oscillatory in angular dependence as 
Piicose) ~ sin[(Z + \)e + f )]/\//sh^ (for W > 1). 

Now our task is to extract useful pieces of information from Eq. (fT2|) with the high-energy 
elastic S'-matrix of Eq. (ITS]) or Eq. (HM . 



The inputs are the isospin-current couphng of p and the exchange degeneracy of p and /2. 
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FIG. 3: Energy dependence of ao(s) in the complex plane (Argand diagram). 

IV. FINAL STATE INTERACTION OF TWO COUPLED CHANNELS 

Let us first count how many dynamical quantities are involved in the most general (n x n) 
FSI relation, Eq. (fT2l) . The unitary and symmetric S matrix contains |n(n+l) independent 
parameters; n eigenphase shifts and |n(n — 1) rotation angles of O. To solve for Aj's of 
observable channels in Eq. f|T2|) . therefore, we must feed the ^n{n + l) dynamical parameters 
of strong interaction. This is not sufficient to determine Ai uniquely. Although the FSI 
relation Eq. f|T2l) may look as if it introduced 2n constraints through the real and imaginary 
parts, only a half of them, namely n of them are actually independent. We must provide 
the relative magnitudes of Aa or A^ as an additional input in order to determine the FSI 
phases uniquely. The magnitude of a decay amplitude is determined primarily by weak 
interaction, i.e., the property of decay operators. Knowledge of strong interaction alone can 
never determine multichannel FSI phases. We need to know interplay of strong and weak 
interactions. 

Nobody is capable of tackling this problem for a general value n. We will therefore be 
content with studying the FSI relation first in the simple manageable case of n = 2 and then 
searching a sensible approximation in more complicated and realistic cases. 

Although the two-channel problem is the next to the simplest, there has been no serious 
attempt to study this case in the past, probably with a good reason as we see below. 
Although it may not look much relevant to the B decay of the real world, we have a chance 
to see through general characteristics of coupled channel effects. For instance, how is the 
FSI phase of tttt channel affected by the pp channel ? If one of the charm- ant icharm channels 
such as D^*'^D^*^ strongly couples to the tttt channel, how does this channel affect the FSI 
of the TTTT channel ?^ While simple-minded perturbative calculations have been undertaken 
in the past, we would like to study these questions systematically with the two-channel toy 



® In the case of a single channel, the relation Ai = e^^'M* gives a constraint only on the phase of Ai, not 
its magnitude. In the case of n channels, something similar happens: The phases of Aa for eigenchannels 
are determined when Sij are completely specified, but their magnitudes \Aa\ are not. 

^ For some dynamical reason the branching fraction to pp is an order of magnitude larger than that to tttt. 
The D* D channel has a huge branching fraction because of the robust b ^ c transition. 
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model that incorporates unitarity. 

We can write the general T-invariant S'-matrix of 2x2 with three parameter (|n(n+l) = 3 
for n = 2) in the form of 

S = 



cos 6 — sin 6 
sin 6 cos 6 




cos sin 6 
— sin 6 cos 



cese 



(20) 



where cos 6' and sin^ are abbreviated as ce and se in the second line. Substituting this 
^-matrix in the FSI relation Eq. f|T2l) . we obtain the constraints on the real and imaginary 
parts or the magnitudes and phases of the decay amplitudes defined by^° 



A, 



(J = l,2) 



(21) 



where the phases Ai 2 are the FSI phases (the strong phases) of channel 1 and 2. We have in 

mind j = 1 for tttt and J = 2 for either pp or D^*'^D^*^ of / = 0. The constraining equation 
of Eq. f|T2|) can be written out for the real and imaginary parts as 



hi 
\b2j 



( ReS \mS 
\lmS -ReS 



hi 
\b2j 



(22) 



where ReS" and ImS" are the 2x2 matrices. To be explicit, 

/ ai \ f 025^4 + C2S24 {c2Si - C2S2)c9Se 
(C2<5i — C2S2)cgSg 025^3^ + €262^ 



02 

bi 
\b2j 



S25,cl + S252SI [S25i - S252)CeSe 
V (S2<5i - S2S2)ceSe S25^Sq + S2S2CI 



S25x4 + -5252 -^e 
(S25i - S252)ceSe 
-ic2S,cj + C252SI) 



■(C25i - C252)ceSe 



[S25i - S252)CeSe, 

S2S1S0 + S2S2Cg 
-(C2<5i - C252)ceS0, 
-{C25isl + C2S2CI) J 



0.2 

bi 
\b2j 
(23) 

As we have pointed out above, Eq. fl23l) contains only two independent constraints, not 
four. Indeed, one can show that two eigenvalues of the 4x4 matrix in the right-hand side 
are unity and generate no constraint. The two remaining eigenvalues are —1 and generate 
constraints. 

We now feed our dynamical input Eq. f[T7|) in the slightly different notation; 



S 
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1-e, (|5n|<l) 



(24) 



where e is real but not necessarily very small in magnitude. For our discussion later we 
choose e to be ~ 0.5, as suggested by the Pomeron dominance in elastic scattering. It is 
easy to include the nonleading Regge contributions and relax the condition e* = e. Fixing 
e amounts to setting two parameters in S'-matrix of Eq. fl2U]) so that we are left with one 
parameter out of three. When we fix the Sn component as in Eq. UM . it is more convenient 
to parametrize the S-matrix in the form 



S 



iV2e 



e e 



2ix 



(25) 



Recall that the arbitrary unphysical phases of states have been fixed up to an overall ± sign by the 
symmetry condition Sij = Sji on the S'-matrix. 
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where the angle x is the remaining single parameter of the 5'-matrix. It is related to two 
eigenphase shifts 81^2 defined in Eq. (l20l) by 

. sin^ 25i - sin^ 282 

sm2(di — 02) 



Note that once Su is given, magnitude of the channel coupling 15*121 = -\/2e — in Eq. (j25 
is fixed by unitarity and no longer a free parameter. 



Case of 1^22 - 1| > I'Sii 



Let us first study the case where the partial- wave amplitude ao'^(s) of rescattering is 
stronger in the second channel than in the first channel, namely, i.e., \S22 — 1| > \Sii — 1|. 
This may serve as a toy-model of vrvr and D*D The D*D scattering of / = is presum- 
ably strong because of the near-threshold enhancement and/or broad excited charmonium 
resonances. For illustration we consider the extreme case that the rescattering in the second 
channel is maximally strong relative to that in the first channel. This is realized by choosing 
e^*^ = — l(x = f ) in Eq. fl25|) . The symmetric unitary ^-matrix takes the form of 
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In terms of partial-wave amplitudes, this gives a'^'^^s) = and i{l — |e) ^ ^ie for the first 
and second channels, respectively, with |e ~ 0.25. By substituting this S'-matrix in the FSI 
relations, we obtain (see Eq. fl2T]) ) 

(28) 




The real and the imaginary parts of Ai are still independent of each other. The phases Ai^2 
cannot be determined uniquely even after the S'-matrix is fully specified. But the phase Ai 
is related to A2 by 

A ^ A 

tanAi = tanA2 

2-e 

~ — X tanA2. (29) 

3 

Even when x of S22 = \S22\g'^''^ is equal to |7r (normally called as "resonant"), A2 is not 
necessarily equal to ^tt. Although this may look puzzling at the first sight, it is not. To 
determine Ai 2 uniquely, we need to feed one more piece of information. For instance, if 
the value of the ratio IA2/A1I is supplied, we can determine Ai and A2 individually. By 
eliminating A2 from Eq. fl29l) we obtain the relation that determines Ai in terms of |yl2/v4i| 
and e: 

M2\ 



sin^ Ai ^ 



4(1 -e) 



|2 



(2 



l^il 



(30) 
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The D*D branching is larger than that of DD. The D*D and DD channels cannot make 
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The ratio contains information of weak interactions. In multichannel decay, weak 

interaction plays a very important role in determining the FSI phases. 

There is one shortcoming of this two-channel toy model: As one sees in Eq. (!30!) . the 
ratio 1^2/^1! must lie in the range of 

(-) 

for this model to be applicable. At the lower boundary of IA2/A1I, it happens that Ai = 
A2 = 0, while Ai = — A2 = ±90° at the upper boundary of |yl2/Ai|. When |yl2/Ai| is in 
between, both Ai and A2 take nonzero values even though all elements of the S'-matrix are 
real. This is an important point to be emphasized. The phase of Ai can arise from the 
process 5 — 2 — 1 through the intermediate state 2. Some may wonder why the ratio 
1^42/^1! is constrained in the two-channel toy model. With channel coupling present, one 
channel feeds the other by FSI to the direction to equalize magnitudes of \Ai\ and \A2\. 
The FSI not only generates phases for Ai and A2 but also alters their magnitudes. Highly 
asymmetric \Ai\ and IA2I are incompatible with the FSI connecting the two channels unless 
e — s> 0, i.e., \Si2\ 0. 



B. Case of ^22 ^ 

When channel 2 is pp, we expect that the elastic pp scattering is asymptotic at y/s ^ 5GeV 
and very similar to the elastic vrvr scattering: 

522 = ^11 = 1-6, (e* = e) (32) 

In this case unitarity and symmetry require that the off-diagonal element 5*12 should be 
purely imaginary (x = in Eq. (1251) ): 



This is the case that was studied by Donoghue et al^. The FSI relation Eq. f|T2l) leads us 
to 

Ai = ai + ibi 
A2 = ^{b,+^^^a, 

tanAitanAa = — ^ ^ -. (34) 

Here again the FSI phases are uniquely determined only after the ratio |742|/|Ai| is given. 
The phase Ai is expressed in terms of |742/v4i| and e by the same relation as Eq. flSU]) . The 
ratio |742/v4i| lies in the same range as Eq. flHTl) . In contrast to the case of the maximum 
l'S'22 — ll, the amplitude A2 and Ai are 90° out of phase, one real and the other purely 
imaginary, at the lower and upper boundaries of the range for IA2/A1I. 
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C. General two-channels 



Once we have explored the two extreme cases above, it is not difficult to find the general 
solution for an arbitrary value of x in Eq. fl25l) . By rewriting the FSI relation in A'2 = e~'^^A2, 
we can reduce it to the case of S22 = Su above. We obtain the solution for general x as 



Ai 
A2 

tan Ai tan( A2 — x) 



ai + ibi 



2 - e 



(35) 



The first and second lines of Eq. (l35l) reduce to those of the previous cases; 1 1 — »S'22 1 — max 
and 5*11 = 5*22 as x l^r and x 0, respectively. The third line relates the FSI phases 
Ai and A2 to each other with parameter e and x- The ratio IA2/A1I is related to Ai and e 
exactly in the same way as in Eq. (!30|) . Consequently | ^2/^1 1 is also restricted to the same 
range, Eq. (1311) . 

We can actually find the exact solution even when we include the very small imaginary 
part of 1 — e in Su due to the low-ranking Regge trajectories in elastic scattering. In this 
most general case it is convenient to write the S'-matrix in the form 



S 



e e 



g2gj(xi+X2) 



(36) 



Unitarity and symmetry fixes the phase of 5*12 as shown above once those of Su and ^22 
are given. For this reason one earlier literature 15|] assigned the angles Xi and X2 to the 
"initial-state" and "final-state" interactions of scattering, hinting that the phases of decay 
amplitudes Ai and A2 acquire Xi and X2, respectively, by FSI. Unfortunately this interpre- 
tation was wrong. As we see below, the phases of Ai^2 are xi,2 P^us additional contributions 
that depend on mixing of the channels and weak interaction. It should also be pointed out 
that the simple phase relation like arg Sij = arg (—iySuSjj holds only for the S-matrix of 
2x2, not of more than two channels. 

To solve the FSI relation with Eq. (1361) we factor out the "elastic phases" xi and X2 of 
the diagonal S'-matrix elements Su {not of the partial- wave amplitudes a/(s)) from the decay 
amplitudes by introducing A'^ by A'^ = e~'^^^Ai {i = 1,2). Then the FSI relation A' = SA * 
reduces to the form identical to Eq. (133|1 . Therefore the solution for 2 can be immediately 
written as Eq. 



A2 



bi 



-ai 



tan(Ai - xi) tan(A2 - X2) 



2-e' 



(37) 



where ai and bi are real The first and the second line of Eq. (!37|) requires |y42/Ai| to remain 
in the same range as in the previous two special cases, (cf Eq. (13T]) ). The relation of Eq. 
(130|) is trivially modified as 



sin^(Ai - xi) 



4(1 



(38) 
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This relation is the bottom hne of the general two-channel toy model: The total FSI phase Ai 
of Ai is sum of the rescattering phase xi of y/Sn, not of the elastic partial-wave amplitude 
{ai = — 1)), plus the the phase due to rescattering through the second channel. It 
cannot be over-emphasized that in the presence of inelasticity the phase of y/Sn is very 
different from the phase of the partial- wave amplitude. For Pomeron-dominated scattering, 
for instance, argao(s) = ^vr since ai{s) is purely imaginary for all /, but argv^STi = or 
±7r since is real and positive for < Imoo < 0.5. The phases of and a;(s) would 
be equal only in the elastic limit where Su = e^^^^ and a/(s) = (l/2i)(S'ii — 1) = e^^^ sin(5i. 
It makes no sense whatsoever even as an approximation to equate the FSI phase with the 
phase of the elastic partial-wave amplitude a; in B decay. 

It is worth mentioning here that the solutions of the two-channel problem, Eqs. (137|) 
and (138|) . apply to the Q~ decay into AK~ and (S7r)7=i/2- The AK~ and En yields add up 
to over 99% of the observed nonleptonic final states. In this case the lower ranking Regge 
trajectories contribute more to the diagonal S'-matrix elements than in the light hadrons 
channels of B decay. 

While the two-channel toy model casts light on many important issues, it has one unde- 
sirable feature that the ratio |y42/Ai| is restricted within the rather narrow range set by Eq. 
( 13T|) . Numerically, 

0.58 < IA2M1I < 1.73. (39) 

This constraint limits applicability of the two-channel toy-model to B decay modes. We 
must extend to more than two channels to study B decay. However, as the coupled channels 
increases, the number of dynamical unknowns quickly increases in the FSI relation. In 
order to keep our problem manageable, we must introduce some approximation that keeps 
mathematical complexity under control. 

V. TRUNCATED MULTICHANNEL PROBLEM 

Going back to our fundamental equation Eq. (fT2!) . we consider the situation where one 
inelastic channel makes a dominant feed back to the elastic channel (channel 1) and all other 
inelastic channels are either unimportant individually or largely cancel among them. We 
are specifically interested in the case, 

l^ail < \Sn\, but l^al > \Ai\ (40) 

such that 

\SuA;\=Oi\SuAl\), lY^S.jA^l-^lS.^All (41) 

i>3 

In this case we can truncate the sum over the inelastic channels at j = 2; 

A,= Si,A* + J2 Si,A* ^ SuAl + SuA;. (42) 

J=l,2 i>3 

The FSI relation for the channel 2 reads 

A2 = S21AI + S22A; + S2jA*. (43) 

i>3 

Since unlike the off-diagonal S21 the diagonal S-matrix element ^22 = 1 + 2ial^{s) contains 
the term unity, we expect that 5*22 is 0(1) or a substantial fraction of it unless the term 1 is 



13 



cancelled accidentally by 2209^ (s) with high accuracy. In comparison, 5*12 represents a small 
leakage into a dominant inelastic channel in the present case. Therefore 

\S2iAys22Ai\ = \S2i/S22\ X (44) 

where the both factors in the right-hand side are small. In other words, the effect of the 
single elastic channel back on a robust inelastic channel j is negligible when the coupling 
5*^1 between them is feeble. Therefore, magnitude of the first term 5*21^^ in the right-hand 
side of Eq. fllSl) is much smaller than that of 15*22^21 by the assumptions made in Eqs. fHD]) 
and dHJ. Therefore, we may drop the first term in Eq. fH51) . Then no information of the 
channel 1 is needed to solve Eq. fH51) for A2. Therefore we solve only Eq. fH2]) and obtain a 
relation between Ai and A2. Solving Eq. (1431) for A2 may be hard. In the numerical exercise 
later we do not attempt to compute for A2 theoretically in terms of other inelastic channels, 
but resort to experiment for information of A2. 

The S'-matrix now need not satisfy unitarity in the subsector of channel 1 and 2. It can 
be written in general as 

S = ( ,!^,^Jf.Z) '"'^t-^r"^ ) > (0 < - < < A < v/T^), (45) 

where two real parameters n and A have been introduced to describe inelasticity of scattering. 
When the channel 2 is also a two-body light-hadron channel, the value of A is ~ 1 — e and 
X2 — Xi- In decay the branching fractions to two-body light-mesons are much smaller 
than those to charmed meson pairs by the property of weak interaction. The decay B Ktt 
through D*D* is a typical example since \k\ is much smaller than 1 — e but |y42| is much 
larger than \ Ai\. It has been speculated that the presence of the D*D* channel may generate 
a large FSI phase for KvrfTij. We will examine this possibility later. 

The FSI relation can be solved for 2 even in the presence of the additional parameter 
K by the rephasing technique that we have used earlier. When we express Ai and A2 in 
terms of 02 and 62 instead of ai and bi, the solution of Eq. f|T2l) with Eq. fH5l) is: 

Ai = ne'^^ Q [-02 sin(x2 + Xk) + &2 cos(x2 + X«)] + ^ZT^ ['^2 cos(x2 + Xk) + &2 sin(x2 + Xk)]^ 
A2 = a2 + ib2, (46) 

As we have pointed out, the parameter A is not needed to express the relation between Ai 
and A2 in the truncated approximation. Although Xk enters Ai, we can express the phase 
Ai without Xk by using experimental knowledge of |y42/v4i|. Rewriting Eq. P6l) with the 
magnitudes and phases, we obtain a simple generalization of the previous relation, 

-^(^.-^O^j^(k^^-.^). (47) 



The right-hand side of Eq. fl47j) gives the contribution of the channel coupling that is to be 
added to the elastic contribution xi- The ratio IA2/A1I is now bounded as 

e 1^2! 2-e , 

If K is small, that is, if the leakage into the channel 2 is small, large values can be accommo- 
dated for 1^2/^1 1- Therefore the truncated model is applicable to more general situations 
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than the two-channel toy model that we have discussed. The ratio |y42/Ai| contains informa- 
tion of weak interaction. It is amusing to see in Eq. (H7!) that the strong phase Ai coincides 
with the small "elastic phase" xi of V Su when A2 takes the smallest value, |y42/Ai| = e/n, 
in the allowed range of Eq. (HSj) . The channel coupling effect |Ai — xi| is the strongest when 
/All takes the upper limit (2 — e)//t in Eq. (HHll . 

If we proceed further and include two prominent inelastic channels, the FSI relation for 
Ai is: 

~ SnAl + Si2A; + Si^A;. (49) 

If we continue along this line and incorporate more inelastic channels, the FSI equation for 
channel 1 turns into 

ai + ibi = Su{ai + i&i)* + SijA*. (50) 

The first term in the right-hand side can be viewed as counteraction of elastic rescattering. It 
affects not only on the phase of Ai but also causes long-distance enhancement or suppression 
on magnitude of Ai depending on the force in the elastic channel. With Su = 1 — e, Eq. 
(I3U]) reveals an interesting general feature of the multichannel FSI. Separating the real and 
imaginary parts of Eq. fl5Up . we have in the case of real Su 

ai = -ReY^SijA*, 

bi = ^T^ImE^i.^l- (51) 

^ e j>2 

Eq. flSTj) shows that elastic rescattering 6*11 enhances the inelastic rescattering effect by 1/e 
(~ 2) for the real part ai of Ai and and suppresses it for the imaginary part bi by 1/(2 — e) 
(^ |). This characteristic depends only on ^n, which we belive we know fairly accurately. 

Despite its simplicity Eq. fl5T]) contains useful information. For instance, if the transition 
to and from channel 1 can be described by the Born terms of t- and w-channel exchanges, 
the off-diagonal partial-wave amplitudes are real so that the off-diagonal S'-matrix Sij {j > 
2) = 2iag^"'^ are purely imaginary. Therefore the real decay amplitudes Aj {j > 2) of 
the inelastic channels contribute to the imaginary part bi of channel 1 in this case. Even 
if there is no resonance in the process, the phase Ai can be very large in this way. If 
furthermore the "inelastic" decay amplitudes Aj (j > 2) happen to be all real, the phase 
Ai would be ±90°. This is not surprising: In the language of dispersion theory the on- 
shell inelastic intermediate states generate an absorptive part for Ai that turns out to be 
purely imaginary in this situation. Our Eq. fl5ip involves one dynamical input: The elastic 
scattering amplitudes of light mesons are almost purely imaginary (Pomeron-dominated). 

VI. NUMERICAL EXERCISE 

Some numerical exercise is called for to show relevance of our endeavor to the B decay 
of the real world. Contrary to the initial optimism that had prevailed before B physics 
experiment started, analysis of experiment seems to suggest that the FSI phases of some two- 
body decay amplitudes appear to be much larger than what we expected in the original short- 
distance picture [23]. The word Kit puzzle has been coined for the unexpectedly large tree- 
contribution and/or FSI phases in Kn modes. As the perturbative technique has become 
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more sophisticated, people have come to agree that emission and absorption of soft and 



colhnear quarks and gluons plays an important role in many decay modes [2 ll. |22| . Such 
soft constituents contribute to the FSI phases involving long-distance physics. While one 
can parametrizes such contributions in the soft-collinear theory, one cannot evaluate them 
numerically in perturbative argument. Our S"- matrix approach also has its own drawback: 
While elastic scattering of two light hadrons at energy tub has been well understood, we 
know less about their inelastic scattering. Nonetheless we would like to show here that our 
method may be useful in some of B decay modes. 

It is believed that the decay B Kn occurs primarily with the penguin interaction 
~ [(bs){qq) + h.c]. In the penguin process the coupling to channels such as K*p need to be 
studied as a source of the strong phase of the Kir amplitude. However, it has been argued 
that Kn can be produced indirectly with the tree interaction ~ [(bc)(cs) + h.c] as well 
through the charmed-meson-pair states such as DDg and D*^*!^]. The Kn amplitude of 
this process acquires a strong phase different from the direct penguin amplitude. What can 
our approach say about this problem ? 

The branching fractions have been measured for the following two-body channels that 
couple to Kt^ in B^ decay |23]: 



B{K+n-) 




(1.88 ±0.07) X 10 






(1.15 ±0.10) X 10 


B{K7]) 


< 


2.0 X 10~^ 


B{K7]') 




(6.5 ±0.4) X 10"^ 


B{K*'^(P) 




(0.95 ±0.08) X 10 


B{K*+p-) 


< 


1.20 X 10"^ 


B{K*°u) 


< 


4.2 X 10"^ 


B{D-Dt) 




(6.5 ±1.3) X 10"^ 


B{D*-D*/) 




(1.77 ±0.14) X 10 



(52) 
(53) 



-2 



(54) 



All of them can make = 0~^, the spin-parity of Kit. The states Ktt and K*p can be 
either in J = | or / = | while all other modes are only in J = i. 

From the Regge phenomenology on elastic Kn scattering, we have already estimated 
e for the Pomeron contribution in Section III. The p and /2 trajectories with exchange 
degeneracy allow us to estimate the small imaginary part of 5*11. By adding the p and /2 
Regge contributions, we obtain numerically 

r 0.39 xe°o^^ (/=i), 
"^11 - \ 0.46 X e-°-i°^ {1 = 1). ^ ^ 

The p/f2 contributions to the scattering amplitudes T are dominantly imaginary for I = \ 
and real for J = |, as we expect from s to t-channel duality. Therefore they generate a 
larger imaginary part for 5* = 1 ± 2iT in the I = \ channel. This is very different from our 
intuitive picture in the single-channel case. Though we do not attach errors to our estimate, 
errors as large as factor two are possible for the phases in Eq. (1551) . 

Estimate of S12 is less reliable than owing to very indirect experimental information 
and to larger theoretical uncertainties. We eliminate the K*(l) channel from our consideration 
since the Regge residue or the coupling of to vrp is highly suppressed ( "OZI-forbidden"). 
In contrast K*p can couple to Kn without such suppression. While a loose upper bound 
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has been set on K*p experimentally (see Eq. we believe that the K*p channel is 

more important than the channel. We may use the Regge phenomenology to estimate 
5*12 for the K*p state of longitudinal polarizations in the final state. For the scattering 
Kvr — » K*p, the leading Regge poles are uj and a2 which are exchange degenerate. The 
couplings of vrpcj and Tipa2 are known on the mass shells of uo and a2, respectively, from 
low-energy spectroscopy. The corresponding KK* couplings of uj and 02 are obtained by an 
SU(3) rotation. But we need to extrapolate them to the off-shell p and 02 to relate them 
to the Regge residues. This extrapolation is a major source of uncertainty. If we ignore the 
extrapolation, they are at the same level in magnitude as the nonleading Regge contributions 
in Sii] 

1 -0.07 + 0.02^ (/ = |), 
^\ -0.05^ (/ = |). ^^^> 

Using the ratio |yl2/Ai| computed with the measured branching fraction and the upper 
bound listed in Eq. (15^ . we reach the crude estimate, 

|^iTM;|<0.1x|5nAt|, (57) 

Despite large uncertainty of these numbers we may conclude with Eq. fl56|) that the inelastic 
term 5*12^2 ^^e K*p channel in the Ktt mode is not significant relative to the elastic term 
SiiA\. It is certainly not a major source of the strong phase for the Kt^ amplitude in 
the standard penguin decay (not through cc). The same line of estimate suggests that the 
K*uo state is neither important to the FSI phase of Kit. Some may wonder about one-pion 
exchange in Ktt ^ K* p. The Reggeized pion exchange amplitude for Kt: — > K*p is down 
by another power of s°'^ relative to those of uo and 02 exchanges. The small denominator 
of the pion propagator does not enhance the amplitude near forward direction because the 
Lorentz structure requires the amplitude near the pion pole to be proportional to 



m2-2(p^-Pp) 



The both factors in the numerator vanish up to 0(4m2/m^, 4m|^*/m^) for longitudinally 
polarized K* and p in the forward direction (pp || p^). They eliminate a forward peak 
from the pion pole. The Regge theory predicts that the rest of the scattering amplitude 
falls sharply (~ e^'^^^-'*''^*) off the forward direction. Therefore the pion exchange can be 
dismissed. Then we feel safe to conclude with Eq. flFTI) that the coupling of Kti to the K*p 
channel is not important in determining the FSI phases for the Kir modes. 

The Krj and Krj' channels can couple to the Kt: channel of / = | through the 02 Regge 
exchange. Although this contribution has the same uib dependence as the p and /2 Regge 
exchanges, the 02 Regge residues with tttjIt]') are most likely smaller than the residue with 
vrp: We can estimate from the 02 decay branching that, after the d-wave phase space factor 
is separated, the on-shell 02 couplings to nr] and tti]' are about factor 20 smaller than those 
to vrp. Therefore neither Kt] nor Kr]' compete with K*p in the final state. Therefore we 
may leave out Krj and Kt]' from our consideration. 



The K* trajectory generates a backward peak is generated, but it is less important than the forward 
peak. We neglect the backward peak contribution to ■ 
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The contributions from nonresonant three-body final states are harder to estimate since 
computation of S12 is next to impossible. There are many nonresonant multiparticle channels 
with relatively minor branching fractions. If rescattering of Kn to multibody channels is 
quasi- diffract ive with no quantum number exchange, the Pomeron can contribute. In such 
scattering the final states consist of two lumps of relatively small invariant masses that carry 
the same flavors as K and vr. They are likely to end up in two-body states of highly excited 
meson states, for instance, i^2(1430)ai. While this is a possibility, none of such modes have 
been positively identified so far in measurement. 

Genuine nonresonant three-body channels are probably not a major source of the FSI 
phases, unless their contributions add up by constructive interference to a large value. In 
fact, it is conceivable that they sum up in random phases into relatively a small numberf^, [^. 
That is one motivation when we have introduced the truncated approximation. Our tentative 
conclusion on the Ktc amplitudes (/ = |, |) of the light-quark penguin decay operators is 
that the FSI phase produced by coupling to the inelastic channels is insignificant. Analysis 
of the Kn amplitudes in search of the weak phases was started more than ten years ago. 
With little knowledge of the strong phases, however, the analysis could be carried out only 
by assuming that the long-distance strong phases be negligible j24|. 

The charmed meson-pair channels are very different. Since they are the CKM dominant 
tree-decay final states, their branching fractions are a few orders of magnitude larger than 
that of the penguin-dominated decay. They can annihilate into Kn. In the quark picture this 
process can be viewed as the on-shell contribution of the cc penguin to Kir. Some call this 
process as "charming penguin" [7|. Among the charmed meson pairs, the most prominent 
decay channel of = is D*D*. Its amplitude can be estimated with Eq. fl54p as 

I^D*-D*+Mi^-(/=i/2)l - 25. (59) 

One important question here is how much of the observed total AKn{i=i/2) is the "charming 
penguin" contribution. Since we expect that long-distance physics enters the on-shell process 
of charmed-meson pairs at energy ttib, it is not easy to evaluate its magnitude. Some argue 



that it can be very large [171. Il8l |. However, a counter argument was made to advocate the 



short- distance argument (l9l|. Theorists have not come to consensus on magnitude of this 
contribution. Therefore we insert one fudge factor r here for this contribution to Ai as 

\AD-Di+/^K7TiI=l/2viacc) \ ^ 25 X i, (r < 1) (60) 

where r, the fraction of the cc contribution, may be as large as a half or even more H. It 
must be settled by theory rather than by experiment. With this fudge factor Eq. P71) turns 
into 

where k' = n/r = |S'i2/tS'ii|. Even if the spill-over of D*D* into Kir is as tiny as one tenth 
of percent {k^ = ISu/Sul"^ ~ 10~^), the D*D* channel may control the FSI phase of the 
Ktt amplitude that comes through cc. To proceed further, we must look into the transition 
Ktt ^ D*Dl. (Fig.4) 

Application of the Regge theory is questionable to the charm-pair channels since the 
charmed meson masses are around 2 GeV and the total energy is a little above 5 GeV. 
Departing from the Regge theory, let us make the Born approximation in t-channel exchange. 
The exchanged mesons are the charmed mesons D, D*, ■ ■ ■. For D-exchange, we know 
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FIG. 4: The dominant inelastic scattering Kit D*Dl. 

the D*Dtt coupling from the decay D* Dtt and can compute the D*DsK coupling by 
an SU(3) rotation of the D*Dtt coupling. The differential cross section rises toward the 
forward direction when the D meson is exchanged. For D* exchange, we can deduce the 
D*D 71 coupling from the D*Dt[ coupling with the heavy quark spin symmetry and rotate it 
into the D*D*K coupling by SU(3). In contrast to the D exchange, the Lorentz structure of 
the vertex ~ £^'^'^^€i,pikP2\ of D* exchange cancels a forward peak that would be otherwise 
generated by the D* propagator. For this reason the D*-exchange is less important. When 
we compute the D contribution in the Feynman diagram with the on-shell couplings, we 
obtain |S'i2| ~ 0.5. But this is obviously a nonsense. The reason is that we have ignored 
the form-factor damping effect of the exchanged off-shell D. For an order-of-magnitude 
estimate we may multiply a factor of ml/mj^ as a form-factor effect where m^, ^ 0.3 GeV 
is the binding scale of the charmed mesons. Then our estimate goes down by more than an 
order of magnitude from |5'i2| ~ 0.5 to |S'i2| ~ 0.014. This latter value is probably closer 
to reality. It is roughly in line with the rule of thumb; in the quark picture a pair creation 
probability of cc is suppressed by about (mg/mc)^ relative to light-quark pair creation of 
qq, where the quark masses are the constituent masses. This rule works roughly for ss and 
cc production in high-energy collision. If we use this rule, we obtain 15*121 — 0.01 from Eq. 
fl56|) with ms/rric — 1/3. Therefore we choose as our best guess 

Si2^*' ~ O.Oli. (62) 

As we have noted earlier, Su is purely imaginary for real ao{s). The number of Eq. fl62|) is 
obviously an order-of-magnitude estimate at best. With e ~ 0.5 and tentatively r ~ 0.5, a 
crude central value of our estimate for n' is 

k'^0.01/[0.5x (1-0.5)] = 0.04. (63) 

We now substitute all these numbers in Eq. (!6T|) of the truncated approximation. We take 
the number of Eq. (163|) as a ballpark figure and sweep the value of /t'^ by a factor two across 
this value to see what FSI phase can be generated for Kn of / = | by the channel coupling 
to the D*D* channel. The result is plotted in Fig. 5. The value of k'|A2/Ai| is constrained 
between 0.5 and 1.5 by Eq. (16T!) and sweeps in the region between two vertical broken lines. 

Within the uncertainty of k', the FSI phase Ai of the Kn amplitude through cc can be 
any value between 21° and 69°. Eq. (I6T1) does not determine the sign of Ai since it does not 
contain full information of A2. Because of the large uncertainty of k', we cannot constrain Ai 
meaningfully at present. Even Ai ~ 90° is not reliably excluded. Keeping the uncertainty 
of our estimate in mind, we should state here only that the channel coupling to D*D* is 
capable of producing a very large strong phase for the charming penguin Ktt amplitude, in 
particular in the case that the "charming penguin amplitude" is a sizable fraction of the 
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FIG. 5: The FSI (strong) phase of the Kir decay channel of / = 2 as the K-k D*D* transition 
ampHtude is varied in magnitude. 

total amplitude. Quantitatively reliable computation of the FSI phase will be possible after 
we have obtained a better theoretical estimate of 5*12 for Kir ^ D*D* as well as magnitude 
of the decay amplitude through cc. Until then we must not set this strong phase to zero 
but leave it as an unknown parameter to be determined by fit to experimental data. If such 
experimental fit clearly requires a large FSI angle for the Kn mode of / = | but not of 
/ = |, we shall be able to assert that the on-shell cc intermediate state plays an important 
role in the decay B —>■ Kit. 

The statement above can be made for many other two-body light-hadron modes. The tttt 
mode couples to pp whose branching fraction is nearly an order of magnitude larger than 
that of TTTT. But perturbative calculation of the vrvr ^ pp transition[16] by the Born diagrams 
without off-shell damping can easily overestimate it. With an estimate of 5*12 along the same 
line as for Kn ^ K*p, the right-hand side of Eq. (1471) comes out to be negative in the case 
of TTTT pp, It implies that the branching fraction of pp is not large enough to affect the FSI 
phase for the vrvr amplitudes of the tree decay ~ (bu){qq) + h.c. (g = light quarks). The tttt 
of / = can be fed also by the charmed hadron channel D*D whose branching fraction is 
two orders of magnitude larger than that of vrvr. For the vrvr amplitude through cc, coupling 
to D*D is the most important source of the strong phase. To obtain the strong phases and 
to compare with experiment, we again need to know about relative importance of the two 
classes of decay; b uud and h ccd for B ^ rni. 

In contrast, the tree-dominated CKM-favored decay modes such as B ^ Dn have no wide 
open inelastic channels. Transition to the D*p channel is as insignificant as the transition 
to Ktc — > K*p in the Kn mode. Since the D*p decay branching fraction is comparable 
to that of Dtt and Su is much smaller than Su (cf Eqs. ( 156|) and (1571)). of its elastic 
scattering, we expect that 15*12^421 ^ l^nA^I and therefore the channel coupling contribution 
is unimportant. It means that the FSI phases are small in these modes and that simple 
short- distance calculation of the phases can produce an answer not far from reality. 

VII. COMMENT AND DISCUSSION 

One of our purposes is to solve the most general two-channel model exactly and to 
clarify the mechanism of generating strong phases in this toy model. We have seen above 
that existence of competing channels completely changes the strong phase from that of 
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the "elastic rescattering phase" in all cases of two-channel 5'-matrix. Unitarity plays an 
important role here. The other purpose is to introduce a feasible approximation scheme 
which may be applicable to cases of special interest in B decay. We have applied this 
method to the decay B Ktt and have made semiquantitative analysis. But its outcome 
is not numerically satisfying because of limitation in available knowledge about off-diagonal 
scattering. We have extended our analysis to the general multichannel case and come to 
one simple interesting observation: Though it may sound odd, inelastic scattering tend to 
enhance the strong phase of the elastic channel most when inelastic amplitudes are real [i.e., 
Sij = iTij = imaginary) rather than imaginary, if the strong phases of the inelastic decay 
amplitudes are small, (cf Eq. flSTl) .) 

It is a big challenge to go beyond the two-channel approximation. As the number of 
channels n increases, the number of parameters in the FSI relation increases as |n(n + 1) 
for strong interaction. In addition we need piece of information from weak interaction. In 
the truncated approximation we have kept only a single dominant one among the inelastic 
channels. Our assumption is that all other inelastic channels are less important or else sum 
up in random phases to become numerically insignificant. If the higher inelastic channels 
do not sum randomly, there must be some underlying dynamical reason for it. In such a 
case some approach orthogonal to ours may have advantage. For instance, it is an approach 
based on quarks and gluons instead of hadrons. 

The major source of the strong (FSI) phase in the multichannel case is the transition 
to inelastic channels. Consequently accurate computation of the FSI phases depends on 
knowledge of the transition S'-matrix at energy between a channel of our interest and 
dominant inelastic channels. That is, we need to know dynamics of long and intermediate 
distances at this energy. The present author is of the opinion that quantitatively we have 
a little better handle on hadron physics numerically than on the soft and coUinear quarks 
and gluons in this territory. But opinions probably divide among physicists of different 
generations. 
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APPENDIX A: REAL ORTHOGONALITY OF TRANSFORMATION 

Expansion of the observable state in the 5- matrix eigenstates |a°"*) is defined by 



At this stage the matrix O is assumed to be only unitary, not necessarily orthogonal. Make 
time-reversal on the scattering amplitude from eigenstate \a) to observable state \j) in our 
phase convention of states under time reversal: 



\r') = T.o,a\an- 



(Al) 




(A2) 
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Operating 5'''^ on Eq. (lAip . we obtain 

ir) = EO,a|a-). (A3) 

a 

Substitution of Eqs. (jM]) and (jMD in Eq. ([M]) gives us 

^0.^^2.5.^ (A4) 

which proves that O is orthogonah 

O*, = 0,a. (A5) 
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